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Abstract 

After briefly reviewing the definitions of classical probability densities for position, Pcl{x), and 
for momentum, Pcl{p), we present several examples of classical mechanical potential systems, 
mostly variations on such familiar cases as the infinite well and the uniformly accelerated particle 
for which the classical distributions can be easily derived and visualized. We focus especially on 
a simple potential which interpolates between the symmetric linear potential, V{x) = F\x\^ and 
the infinite well, which can illustrate, in a mathematically straightforward way, how the divergent, 
(5-function classical probability density for momentum for the infinite well can be easily seen to 
arise. Such examples can help students understand the quantum mechanical momentum-space 
wavefunction (and its corresponding probability density) in much the same way that other semi- 
classical techniques, such as the WKB approximation, can be used to visualize position-space 
wavefunctions. 



1. Introduction 



Since the introduction of quantum mechanics early in the last century, practitioners and 
students of the subject alike have used semi-classical methods to calculate, understand, 
and visualize properties of position-space solutions of the time-independent Schrodinger 
equation. WKB-type solutions for example, provide an attractive and mathematically 
understandable approximation which makes manifest the typical correlations between the 
magnitude of a stationary state wavefunction and the corresponding 'wiggliness'. At a 
somewhat less sophisticated level, many textbooks , make direct comparisons between 

the quantum mechanical probability density, defined by Pqm{x) = \ip{x)\'^, with intuitively 
derived classical counterparts, Pcl{x), which are often motivated by very simple ^how much 
time is spent in a given interval^ arguments. 

In comparison, there are far fewer discussions in the pedagogical literature or in stan- 
dardly used textbooks giving examples of the corresponding semi-classical connections to 
momentum-space solutions, (pip). Solutions to the Schrodinger equation obtained directly 
from a momentum-space formulation (such as for the uniformly accelerated particle) are 
sometimes presented, but most often the information encoded in (f){p) is simply noted to be 
related to that in the position-space wavefunction by a rather formal operation, namely the 
Fourier transform, via 

(j){p) = ^= ^{x)e'P''/''dx. (1) 

Correlating the information clearly evident in the magnitude/wiggliness of a position-space, 
ipi^), with the similar information encoded in is a non-trivial task, even for experienced 
scientists, and students first approaching the subject often can find no obvious connections 
between the form of the quantum mechanical probability density for momentum, 
and their intuition about classical mechanics. A recent study of the progression of student 
understanding of quantum mechanical concepts, starting with introductory modern physics 
courses, through advanced undergraduate quantum mechanics classes, to first year graduate 
quantum theory courses, has provided evidence that while student understanding of the 
form of position-space wavefunctions does increase through the undergraduate curriculum, 
their ability to extend these notions to momentum-space ideas lags far behind. 

Finally, since the results of experiments from a wide variety of areas in physics are in- 
creasingly presented in the form of momentum distributions jsj, j^, a background in un- 



derstanding a momentum-space approach to classical and quantum mechanics is becoming 
more important for students to acquire. 

Motivated by all of these factors, in this note we will discuss several simple model systems 
for which one can rather easily calculate, analyze, and visualize the classical probability 
densities for both position and momentum, and to also compare them directly to their 
quantum mechanical equivalents to visualize the correspondence principle limit in both x- 
and p-space. We first briefly review, in Sec. 2, the basic definitions of classical probabihty 
densities for one-dimensional bound state problems and apply these methods to the simple 
problem of the 'bouncer', a particle subject to a uniform gravitational force, but with a 
rigid, infinite wall constraint due to a horizontal surface. As our main example, in Sec. 3, we 
provide an attractive interpolating case which shows how both the classical and quantum 
mechanical descriptions for the symmetric linear well (defined by V{x) — F\x\ and hence 
related to the 'bouncer' or uniform acceleration case) and the standard infinite well potentials 
can be smoothly connected to one another, and how the classical momentum distributions 
can be easily extracted, we also show, in some detail, how the magnitude/ wiggliness of the 
position space quantum wavefunction is directly correlated to the form of its momentum- 
space (j){p) transform. This last case also provides an intuitively clear and mathematically 
appropriate description of how the limiting case of the divergent, 5-function momentum- 
space classical probability densities can arise for the infinite well case. 

2. Review of classical probability densities for position and 
momentum: the 'bouncer' as an example 

In order to review the calculation and visualization of classical probability densities for 
position, Pcl{x), and momentum Pciip), we begin by illustrating, in Fig. 1, a rather generic 
potential energy function, V{x) versus x, which can give rise to bound states which can be 
described either classically or quantum mechanically. A fixed value of the energy E is shown 
which then defines the classical turning points, a and b, defined by V{a) = E = V{b): this 
implies that Pcl{x) will be non- vanishing only over the range (a, 6). The corresponding 
range in allowed momentum values are given by conservation of energy using 

p{x) = ±,/2m{E-V{x)) giving = ^2m{E - V^^) (2) 

so that the classical probabihty density, Pcl{p): will be defined over the range (— Pm, +Pm)- 
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The variation of momentum values with position in the potential well is perhaps best shown 
in a phase-space x — p plot, also included in Fig. 1, which simultaneously illustrates the 
allowed ranges in both variables. 

We note that a particle at a given x location will experience both positive {+p{x)) and 
negative {—p{x)) values with equal probability as it moves from left to right and then back 
from right to left. Thus, even if the potential energy function is not symmetric, the classical 
probability density for momentum, Pciip), will necessarily satisfy Pcl{—p) — -Pcl(+p). 
This feature is also echoed in the x—p plot in Fig. 1 where the phase-space curve need not 
be symmetric about any vertical axis, but will necessarily be symmetric about the p = 
line. 

We note from the figure that there may well be more than one location where the particle 
has the same momentum value, given by p = ■\/2m{E — V{x)), such as at x = c and x — d. 
We must then take both contributions to Pcl{p) into account in the formalism which follows. 
One case in which this may not occur is if there is an infinite wall type potential barrier from 
which particles simply 'bounce' and rebound instead of slowing down on one side, coming 
instantaneously to rest, and then reversing direction as shown in the case in Fig. 1. 

The classical motion in any such bound state system will be periodic with period Tcli 
corresponding to a complete set of left to right and back to left (e.g., a to 6 and then back 
to a) motions. For some purposes, we only require the 'half-period' (the a to 6 time) given 
by r = Tcl/ 2 and this quantity can be easily calculated from the potential energy function 
via 



A standard argument leading to the classical probability density for position, Pcl{x).i 
begins by noting that the small probability of finding the particle with position in the range 
(x, X + dx) should be proportional to the amount of time spent in that interval, dt. The 
fraction of time spent in that small interval will be dt/Tci and using the definition of 
probability density, we find that 

Pcl{x) dx = (iProb[a;, x + dx] = — — . (4) 

TcL 

We then recall that (i) the particle will be in the same dx bin on both its 'back' and 'forth' 
trips, so the last term should really be doubled, and (ii) the classical speed, v{x) = \dx/dt\, 




(3) 
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can be used to relate dx directly to dt, giving the familiar result 

dx \ dx J Tv{x) 

where r = is the 'half-period'. We note that Eqn. ensures that this definition will 

give a Pcl{x) which is automatically normalized to unity. 

A similar expression for the small probability that a measurement will find the particle 
with momentum values in the range {p, p + dp) is given by 

Pcl{p) dp = (iProb[p, p + dp] = (6) 

and we need the definition of the classical force, \F\ = \dp/dt\, to be able to write 

dFToh[p,p + dp] dt/TcL 1 

Pcl{p) = ^ = — 1 = ^ ir^. ^| 7) 

dp dp Tcl\F{p)\ 

where we require the magnitude of the force, as a function of p. Depending on the nature of 
the potential, we may have to take two or more contributions into account when evaluating 
Eqn. such as with the x = c,d points in Fig. 1 mentioned above. 

As an example of a calculation and visualization of both Pcl{x) and Pcl{p) for a simple 
potential, we consider the classical 'bouncer', corresponding to a particle bouncing, with- 



out loss of energy, under the influence of gravity, over a rigid horizontal surface 
corresponding potential energy function is given by 



The 



, ^ , l-oo for 2; < 
V{z) ={ . (8) 

mgz for 2; > 

If we imagine a point object with total energy E, the maximum height, if, (which is one 
turning point, in addition to 2 = where it 'bounces'), is given hy E = mgH; the maximum 



momentum magnitude, Pm, is given hy E = p^/2m or pm = v2rnE. The 'half-period', r, is 
simply the time it takes to fall through the distance H, namely r = \J2H/ g and Tcl = 2r. 
The position dependent speed is given v{z) — mgz so that the classical position 

probability density is 



CL 



Tv{z) 



1 



' . (9) 



[2 / m){E — mgz 

which has an (integrable) divergence aX z = H where the particle slows down and reverses 
direction, but is finite at z = where the particle simply 'bounces', reversing direction, but 



with no change in speed. The corresponding classical probability density for momentum will 
be non- vanishing only over the range {—pm, +Pm) and is given by 



Pcl{p) ^ 



TcL\Fip)\ 



1 



2^Wrg 



mg 



1 1 1 , , 

(10) 



2^2m{mgH) 2V2mE 2p„ 



corresponding to a 'fiat' momentum distribution over the allowed range. For this potential, 
due to the infinite wall' at z = 0, the momentum only achieves a particular value once during 
each complete cycle (and not at two separate locations, such as the x = c,d case in Fig. 1), 
so we only have one contribution to |-F(j9)|. (The quantum mechanical solutions, for both 
ip{x) and 0(p), and their comparisons to the classical results, for a very similar problem are 
discussed, in details, in the next section.) 

The expressions in Eqns. (0) and 1111)) for Pcl{x) and Pciip) can be visualized using 
a 'projection of trajectory' technique sl which implements the definitions of the classical 
probability densities in Eqns. (jlj) and © involving 'time spent' arguments. We can use a 
typical solution to the 'bouncer' problem of the form 



z{t) = vqI — gt^ /2 and Vz{t)=VQ — gt (11) 



where Vq = y/2gH is the initial speed, corresponding to an object 'thrown upwards'; the 
motion then repeats itself after each classical period. In Fig. 2, on the respective vertical 
axes {z(t) and Pz(t) = mvz(t)), we indicate equal-sized bins {dz and dp), project them 
horizontally until they intersect the trajectory curves, then 'drop down' to the t axes. The 
small amounts of time, dt, corresponding to each 'bin' in either z or pz are then proportional 
to the probability in that interval and are shown (as a histogram) on the appropriate vertical 
axes. The characteristic peaking near z = H and the finite value at 2; = is evident for 
Pcl{z), while the striking 'flat' momentum probability density is easily visualized using this 
technique. A number {N = 1000) of random 'measurements' of the particle's position and 
momentum over a single classical period are also shown alongside the 'binned' probabilities 
(as dots) and reflect the same phenomena. 

3. The infinite well as a limiting case 

The most familiar of all quantum mechanical bound state problems, the infinite well, 
can also be used an example of classical and quantum mechanical probability densities for 
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momentum. We will also discuss a very useful interpolation between the symmetric linear 
well, defined by the potential V{x) — F\x\ (and hence analogous to the 'bouncer'), and the 
standard infinite well which can be used to show how the position- and momentum-space 
distributions smoothly go from one limit to the other. We first recall that if we define the 
infinite potential with walls at x = ±a, via 

{-|-oo for \x\ > a 
(12) 
for < a 

the classical speed and 'half-period' are trivially related to the total energy via v — y^2E/m 
and T — 2a/v so that the classical probabihty density for position is given 

independent of E and determined solely by the geometry of the well. The energy eigenstates 
for the corresponding quantum mechanical problem can be characterized by parity with, for 
example, the even states being given by 

X = ^ cos ^ — with energies = (14 

^ ^ ' -/a \ a ) " 8ma2 ^ ' 

and similar expressions for the odd states. For large values of n, the quantum mechani- 
cal probability density, Pqm{x) — \'ijjii~\x)\'^ , locally averages to the classical result, since 
(cos2(^)) = 1/2. 

The classical momentum probability density can be easily described in words as consisting 
of equal probabilities of finding the particle with momentum values given by ±po ~ ±-\/2m£', 
which one can write in the form 

Pcl{p) = \ [5{p - Po) + S{p + po)] • (15) 

The corresponding quantum mechanical solutions mimic this 'twin peaks' structure as one 
can see by calculating the momentum-space wave function for the even case, e.g. 

1 /'+°° 



2Trn 

Ap 

2^ 



sin((n — 1/2)71 — ap/h) sin((n — 1 / 2)11 + ap / K) 

{{n -l/2)n - ap/h) ^ {{n - I /2)ti + ap/h) 
sm{{pn - p) / Ap) ^ sm{{pn + p) / Ap) 

{Pn - P) {Pn + P) 



(16) 



where pn = {n — l/2)7rh/a and Ap = h/a. In a classical limit, where we might take ^ — > 
or a oo (or at least much larger than any quantum dimension), we have Ap —>■ and one 
can show that 



lim |0i+)(p)p 

n,/a— >0 



[S{p-Pn) +S(,P + Pn] 



(17) 



as expected. 

To students first encountering such concepts, however, the use of the mathematical for- 
malism of the Dirac 5-function can be less than intuitively obvious, so it would be useful 
to have an appropriate limiting case where students can see, both physically and mathe- 
matically, how the infinite well case is reached. To this end, we consider a variation on the 
symmetric linear well, V{x) = F\x\, namely one with infinite walls included, defined via 



V{x) 



F\x\ = Vo\x\/a for |x| < a 
+00 for \x\ > a 



(18) 



n 



This potential is sometimes called the 'bouncer on a closed court' |9]. In the limit that Vq — *• 
0, this potential approaches the symmetric infinite well. We will especially be interested in 
cases where the energy is greater than Vq (as in Figs. 3, 4, and 5). In that case, because 
of the infinite walls, the particle will never slow down and come to rest as it reverses its 
direction, but rather simply 'bounce' off the infinite wall. In this case, there is not only 



a maximum value of momentum given by = y/2mE, but also a minimum value given 
by p(_) = ^JlmiE — Vq), so that the classical probability density for momenta will only be 
non-vanishing over the intervals (— p+, — P-) and (+P-, +P+)- The 'half-period', r, required 
for the position probability density is given by (so long as > Vq, as we assume) 

doc 



T 



ViX] 



m f Aa 



y/2{E-Vo\x\/a)/m 



(19) 



Vn 



and the classical period is Tcl = 2t. Since the particle achieves the same momentum value 
twice during each period (due to the symmetry of the well) and \F\ = Vo/a, we can write 



Pcl{p) 



2 1 


1 


1 


1 


Tcl\F{p)\ 2 


V2mE - y/2m{E - Vq) 


~ 2 


.Pi+)-Pi-). 



(20) 



over the two allowed regions and vanishing elsewhere. The complete expression for the 
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classical probability density for momentum can therefore be written in the form 

P.art-I^'"^'-'-'^'-"" (21) 
I otherwise 

where Ap = —p{-) and Pcl{p) is clearly normalized appropriately. For fixed values of E 
and a, as we let Vq — > 0, we have Ap — > and we clearly reproduce the (5-function structure 
of the infinite well as the two, isolated rectangular peaks become narrower and higher. 
The classical probability density for position can also be easily derived and we find that 

Pcl{x) = = ^ , , for |a;| < a (22) 

Tv{x) 4a[VE -VE- Vo]^/E - Vo\x\/a 

which has the appropriate limit, namely Pcl{x) = l/2a, when Vq 0. 

In order to compare these classical predictions to the corresponding quantum mechanical 
solutions, we note that the general solutions to the Schrodinger equation for this problem 
can be written in the form 

^(^) = I ClAi{-{x + a)/p) + DLBt{-{x + a)/p) for -a < a; < ^^^^ 
1 CRAi{{x - a)/p) + DRBi{{x - a)/p) for < x < +a 

whe.e Mi,.) a.e ..and„d Ai.y func*. Q and whe.e we have defined 

/ h^a , Ea 

P 



and a= — (24) 

We note that symmetry arguments can be used to relate the L and R coefficients for the 
even or odd eigenstates. The divergent Bi{z) solutions must be used in order to properly 
match the boundary conditions both at the infinite walls and at the origin. Application 
of the boundary conditions at x = and x = a, for example, gives the energy eigenvalue 
condition for the odd states (where ip^O) = 0) as 

and a similar condition for the even eigenstates. The corresponding momentum eigenstates 
can then be obtained by direct Fourier transform from the allowed iplx). 

As examples of how the limiting case of the infinite well can be approached from the 
symmetric linear potential, and of how well the classical and quantum probability densities 
agree (correspondence principle limit), we consider several cases with decreasing values of Vq. 



Specifically, in Figs. 3, 4 and 5, we show examples of solutions (both classical and quantum 
mechanical) for position- and momentum-space probability densities using Eqns. fl21|) and 
(j22|) for the classical versions and the appropriate quantum solutions. In each case, we have 
picked an energy corresponding to an allowed quantum eigenstate with E ^ 10 (in scaled 
units) so that the maximum value of momentum, = \j2mE is virtually identical in all 
three cases. In each case we use fi = 2m = 1 for simplicity, and choose a = 25, while varying 
the values of Vq. We then show results for Vq = 10 (Fig. 3), Vq = 6 (Fig. 4) and Vq = 2 
(Fig. 5) for which the appropriate parameters are shown in Table I. 

We note that as Vq ^ 0, both the classical and quantum position probability densities, 
Pcl{x) and Pqm{x), approach the appropriate 'flat' result for the infinite well: in each case 
we have used the same vertical scale to show this progression. This limiting behavior of 
iIj{x) as Vq — > can be examined more rigorously using 'handbook' properties of the 
Ai{z) and Bi{z) functions. 

For the case of the momentum densities, the approach to the highly peaked case of the 
infinite well is also clear, with decreasingly small values of Ap = — p(„) and increasingly 
large 'spikes': for these cases, in order to see the detailed structure, we have not shown 
the momentum densities with the same vertical scales (note the axis labeling carefully), 
but have kept the horizontal axes identical. Interested students can study the systematic 
variations and correlations between these three cases and may especially notice how the 
spatial variations in the frequency of zero-crossings in the |^/'(x)p are correlated with the 
allowed range in momenta. 

As a final comment, we note that the 'intrinsic' spread of the two prominent individual 
peaks in the momentum space probability densities given by Eqn. ()16|) is of the order of 
Ap = h/a. This means that when the difference Ap = ^(4.) — becomes of this order, 
the classical description will clearly be a bad representation of the quantum results. In the 
numerical studies represented in Figs. 3, 4, and 5, with the specific parameters we have 
used for illustrative purposes, we have Ap = h/a = 1/25 = 0.04 which is still substantially 
smaller than the Ap = p(+) — p(_) values shown in Table I. 

4. Conclusions 

We have examined the structure of classical probability densities for momentum, Pcl{p), 
in several simple variations on familiar, one- dimensional potential problems. As our main 
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example, we have shown how a simple interpolation between the symmetric linear poten- 
tial and the infinite well, provided by Eqn. (fTH|) can illustrate, in an intuitively attractive 
and mathematically appropriate manner, how the divergent ^-function form of the classical 
probability density for momentum can arise in the most familiar of all bound state prob- 
lems, namely the infinite well. Such examples can hopefully provide experience on how to 
understand and visualize the information content of both classical and quantum mechanical 
bound state problems found in a momentum-space description of these systems, a topic 
which is seldom covered in detail in standard undergraduate texts, but which is increasingly 
important in the modern research literature. 
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Tables 





a 


E 






Ap 


10 


25 


10.066 


0.257 


3.173 


2.916 


6 


25 


10.073 


2.108 


3.174 


1.156 


2 


25 


10.105 


2.847 


3.179 


0.332 



Table I. Values of potential well parameters for the 'closed-court' potential of Eqn. p8p used 
in Figs. 3, 4, and 5 respectively. We use h = 2m = 1 in each case. The energy eigenvalues 
are chosen to be as similar as possible, so that the values of are virtually identical in 
all three cases. 
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Figure Captions 



Fig. 1. Plot of a generic potential energy function, V{x) versus x, which supports bound 
state motion. A constant value of the total energy, E, defines the classical turning 
points a, b. The corresponding phase-space x —p plot for the system is shown directly 
below. The extremal values of momentum, ±Pm = ±a/-E — Vmin, correspond to the 
position of the minimum of the potential energy function. The x — p phase-space 
plot is automatically symmetric about the p = axis due to the equal probabilities 
of finding the particle moving through the same dx bin, to the left or right, at two 
different times during its classical period. 

Fig. 2. Illustration of the 'projection of trajectory' technique for visualizing the classical prob- 
ability densities for position (top) and momentum (bottom) for the classical 'bouncer'. 
Also shown are = 1000 measurements of the particles position (z) and momentum 
(pz) coordinate, which are consistent with the binned probabilities shown on the ver- 
tical z{t) and Pzit) axes. 

Fig. 3. Potential energy function (top), classical {Pcl{x), dashed) and quantum {Pqm{x) = 
\ip{x)\'^, sohd) position probability densities (middle), and classical {Pcl{p), dashed) 
and quantum {Pqm{p) = |0(p)P, solid) momentum probability densities (bottom) for 
the symmetric linear well with infinite wall potential in Eqn. (|T8|) with Vq = 10. The 
values of the other parameters are shown in Table I. 

Fig. 4. Same as Fig. 3, but with Vq = 6. Note that the vertical axis for the Pcl{p), Pqm{p) 
plot is not the same as in Figs. 3 or 5. 

Fig. 5. Same as Fig. 3, but with Vq = 2. Note that the vertical axis for the Pcl{p), Pquip) 
plot is not the same as in Figs. 4 or 5. 
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FIG. 3: 
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FIG. 4: 
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l(l)(p)P VS. p 
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FIG. 5: 
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